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The investigation of the decorrelation efficiency of the HMC algorithm with respect to vacuum topology is 
a prerequisite for trustworthy full QCD simulations, in particular for the computation of topology sensitive 
quantities. We demonstrate that for -^-ratios > 0.69 sufficient tunneling between the topological sectors can be 
achieved, for two flavours of dynamical Wilson fermions close to the scaling region (/3 = 5.6). Our results are 
based on time series of length 5000 trajectories. 



1. Introduction 

Topology is a fundamental feature of contin- 
uum quantum field theories with important bear- 
ings on elementary particle physics issues, as 
there is the proton spin content, the large rf mass 
or the role of instantons for the vacuum structure. 

Stochastic lattice sampling of QCD configura- 
tions must be sufficiently ergodic in the topolog- 
ical sectors of phase space in order to allow for 
the proper computation of topological quantities. 
In particular for full QCD, the tunneling rates in- 
duced by the hybrid Monte Carlo algorithm be- 
tween the topological sectors are expected to de- 
crease dramatically when approaching the chiral 
limit, as known from staggered fermion simula- 
tions For full QCD with dynamical Wil- 
son fermions the actual size of the window 
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accessible by todays simulations is still an open 
question which we are going to address here [Q. 

In the last years, various techniques for the reli- 
able computation of topological continuum quan- 
tities from the discrete lattice data have been put 



forth II . In the present investigations, we will 
just make use of the field theoretical definition 
|(| of the topological charge density^ 

Q{X) = 64^^^^)^^ (2) 

with application of cooling Q| as it is not neces- 
sary to compute the renormalization of Q. 

For our investigations, we exploit the SESAM 
and T%L || samples consisting of three time his- 
tories on lattices of size 16 3 x 32 at k = 0.156, 
0.157, and 0.1575, at ^ ratios of 0.839(4), 
0.755(7), and 0.69(1), respectively and of config- 
urations on 24 3 x 40 lattices, again at (3 = 5.6, 
with k = 0.1575 and 0.158, the latter giving rise 
to^« 0.58(2). The length of the samples is 
5000 trajectories. For k = 0.158, the length is 
3500 trajectories. From these samples, we take 
sets of N — 200 'decorrelated' configurations. 

2 Alternatively, we follow Smit and Vink avoiding cool- 
ing. With G being the quark propagator we compute Q 
via stochastic estimates pi: 
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Q = m k p {Tr (75 G))u. 
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Table 1 



First part: mobility D 2 5- Second part: mobility 
normalized by dividing through square root of the 
volume. 
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2. Results 

In Fig. ([!]) we present the time series of the 
topological charge Q from the five samples. The 
first three canvases are computed from SESAM 
ensembles measured on lattices of size 16 3 x 32, 
the next two are from T^L ensembles from 24 3 x 
40 systems. We analyzed every 25th configura- 
tion. In the first three graphs, the results from the 
gluonic and fermionic measurements are super- 
imposed demonstrating nice agreement. For all 
three SESAM masses, the figure tells us that HMC 
can create sufficient tunneling in the SESAM ^ 
window. 

We quantify the tunneling between the topo- 
logical sectors defining the mobility for tunneling 

as 

1 N 

D d = -Y / \(Q^ + 1 )- d )-Q('- d ))l (3) 

i=i 

Tab. (Q) demonstrates that the mobility is de- 
creasing when going more chiral. 

With the 24 3 x 40 system, at n = 0.1575, we 
can assess the volume dependence. Note that Q 
is an extensive quantity. Still we find the tunnel- 
ing rates sufficient, however, with smaller quark 
mass, at k — 0.158, the mobility decreases. As 
we expect the mobility to scale approximately 
with the square root of the volume, we normalize 
the numbers accordingly, see the second part of 
Tab. ©. 

Analyzing the time history on every second tra- 
jectory (at K = 0.1575) we find very frequent tun- 
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Figure 1. Time series of unrenormalized topolog- 
ical charge. We plot series from 3 SESAM and 2 
T%L ensembles. 



neling, see Fig. (||). On this series, we can com- 
pute the autocorrelation function. The exponen- 
tial autocorrelation time is r exp = 80(10). The 
integrated autocorrelation time turns out to be 
Tint — 54(4), compatible with the decorrelation 
of other observables as found in jl0| . 
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Figure 2. Fine scan of unrenormalized Q for k = 
0.1575 on the 16 3 x 32 lattice. 



Histogramming Fig. (g) we find a smooth Gaus- 
sian, see Fig. (||). This is of course less pro- 
nounced for the small decorrelated samples of size 
200, see Fig. (||), however, the distributions are 
still well peaked at Q = 0. 




3. Summary 

As the main result, we find that HMC tunnels 
well between the vacuum sectors for — > 0.69. 

nip — 

Autocorrelation times are slightly larger than for 
the yet worst case observable, the minimal eigen- 
value of the Wilson fermion matrix. Therefore, 
we are confident to be able to carry out reliable 
computations of hadronic properties related to 
topology, like the proton spin content and the 
topological susceptibility for —-values down to 

0. 69. 
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